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Abstract 

A single porous silicon (PS) layer sensor structure was applied to analyze the kinetics of 
adsorption/desorption processes of water, ethanol and ethanol-water solution using CIE1931 chromaticity 
coordinates based parameters as dominant wavelength k dom , complementary wavelength A com and 
excitation purity (EP). Matrix method was used to obtain reflectance spectra dependences on porosity (P), 
PS layer thickness (L P s), refractive index of analyte (n ana iyt e ). Dependences of A, dom , A com , EP on P, L P s, 
iWiyte were computed on the base of computed reflectance spectra and optimal parameters were 
determined for distinction of model analytes as water and ethanol. Time dependences of A, dom , A com were 
experimentally obtained using CCD-based spectrometer for water, ethanol and ethanol-water solution. 
Advantage of characterisation of whole complex multiple peaks spectra using chromaticity parameters as 
opposed to sensor characterization with parameters of individual interference peaks as peak position was 
demonstrated. 
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Introduction 

Porous materials (PM) as porous silicon (PS), are a perspective and promising materials for sensor 
application [1]. PM characterized with large surface to volume ratio in the order of 500 -1000 m /cm [2,3] 
that lets use PM for sensor application. PM properties as porosity (P), pore morphology, thickness can be 
controlled with etching current density and depends on wafer doping level, composition of the electrolyte 
[4]. The penetration of inorganic and organic molecules into pores that leads to changing of optical [5] 
properties of porous layer. Since the effective refractive index of a porous layer depends on refractive 
index of porous matrix and refractive index of the medium inside the pores. Substitution of air inside the 
pores after analyte infiltrated into pores causes an increase of effective refractive index of layer [6] and 



shift of reflectance spectrum to longer wavelength region. Thus, by analysing the reflectance spectrum, 
one can detect the adsorbed molecules of analyte inside the pores. 
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The reflectance spectrum of thin layer has form of a sequence of interference maxima and minima. 
Variation of peak position Z pea k and full width at half maximum (FWHM)) at molecule adsorption can be 
used as sensor parameter[7, 8]. Usually, peaks of reflectance of thin layer are wide enough that degrades 
the accuracy of sensor measurements. In order to monitor more precisely change of spectrum caused by 
analyte, the sensor structure are formed as stacks of thin PS layers such as Bragg mirror (BM), rugate 
filter, Fabry-Perot filter, luminescent PS microcavity [9-14]. Redshift of interference peak maxima 
increases with increasing of refractive index of analyte and can reaches hundreds on nm in infrared region 
[15]. Measurement methods as reflectometry for interference peak characterization and 
photo luminescence can be applied to the same sensor structure and complement one another [16]. 
Sensitivity of single and multiple PS layers sensor structures depends on properties of measuring optical 
systems and for system with spectral resolution of 0.1 nm able to detect refractive index change about 2 x 
10 4 [17] and concentration of different gases and vapor in ppm range [18-20]. Noisiness of spectrums, 
presence of local maximums and not symmetric form of investigated peaks complicates determination of 
/.peak and FWHM, especially for dynamically changing spectrums during adsorption/desorption processes. 
One can expect that the accuracy of the measurement can be increased if the spectrum as a whole is 
analyzed, for example, using the method CIE 1931 chromaticity color space. Now CIE 1931 chromaticity 
color space method is widely used for characterisation of emission spectrums of light emitted structures 
[21, 22]. Electromagnetic visible spectrum in range from A, m i n =380 nm to /. max =780 nm is linked with 
physiological perceived colors in human color vision via chromaticity coordinates as tristimulus values X, 
Y, Z [23, 24] in CIE 1931 chromaticity color space: 
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where S(Z) - spectral reflectance (or transmittance) of sample; 

I(Z) - spectral power distribution of a illuminant. 

x, y, z - CIE defined color matching functions from the color matching experiments [25-27]. 



The color as characteristic of reflectance spectrum can be determined with the x-y parameters as 
coordinates (3) and can be depicted on a two-dimensional graph called the chromaticity diagram. 

X 


x = 


y - 


X + Y + Z 
Y 

X + Y + Z 


(3a) 

(3b) 



x coordinate 

Fig. 1. CIE 1931 chromaticity diagram and computed sample color chromaticity coordinates x-y 
dependence on refractive index of analyte n ana iyte filling the pores (L PS =500 nm, P=50%). 

The outer curved boundary is the spectral locus and each point of locus represents a pure hue of a 
single wavelength. White illuminant in this case has coordinates x=l/3 and y=l/3 on chromaticity 
diagram. Point with these coordinates is called white point (point marked as o in Fig. 1). In the CIE color 
coordinate space for considered color A with x-y coordinates (Fig. 1) there is strait line between this point 
and the white point W. This line can be extrapolated to intersect the border of the locus. In general there 
can be two points of intersection (points D and F on Fig. 1). The point of intersection that is next to the 
considered color on extraporation curve is the dominant wavelength kd om . Considered color always lies 
between white point and point of intersections that determines Adom- The point of intersection on the 
opposite side of extrapolation line relatively to white point defines the complementary wavelength A, com . 










Line DF on Fig. 1 is extrapolation line. Dominant wavelength Ld om and complementary wavelength A com 
are parameters that describe polychrome light mixtures in terms of the monochromatic light that induce 
an identical perception of hue in the human eye [28]. Another parameter that describes color (spectrum) 
in the CIE color coordinate space on chromaticity diagram is excitation purity (EP). EP is determined as 
ratio of segment length AW between white point W and considered color A to segment length DW 
between white point W and of intersection point D on extrapolation curve: EP=AW/DW. EP 
characterizes the saturation of a particular color. The whole spectrum with any number of spectral 
features can be described with only pair of parameters Ldom and EP or pair X com and EP on chromaticity 
diagram. CIE color spaces are applied for color characterization of nanoparticle systems [29], food 
freshness via color change [22], PH monitoring [30], live cell characterization [31]. 

In present work authors analyze kinetics of adsorption and desorption by the measuring of 
reflectivity in single PS layer structure and express received dependences using CIE 1931 color space 
based parameters as /.d« m , A com , EP. Ethanol, water and ethanol-water solution were chosen as matter for 
interest. Dependences of Ld om , A com , EP on L P s, P were computed for wide range of refractive indexes of 
analyte n ana iyt e inside pores. 


Experimental 

Samples and measurements of reflectance spectra 

As subject of investigation the optical sensor based on one porous silicon layer was chosen. PS is 
fabricated using method of electrochemical etching of Si wafer in mixture of aqueous HF solutions and 
C 2 H 5 OH in teflon cell [1,4,7,14]. PS layers with area S=0.79 cm' were formed using method of 
electrochemical etching of p ++ -Si wafer (0.001 Dcm) in HF(40%):ethanol solution in teflon cell. Thermo 
Fisher micro pipette dispenser was used to deposit liquids analytes on the samples. After each 
measurement completion initial state of samples was restored using washing of samples in ethanol and 
drying with nitrogen flow. Control measurements of spectra after washing and drying depict the 
reproducibility of the initial spectra. All spectral measurements were performed with Thorlab CCD 
spectrometer CCS200 at 100 spectrum/s. 

Position of peak maxima X k {n anahte ) shifts to long wavelength region for pores fdled with 
analyte (n analyte ) relatively to position of peak maxima X k (n air ) of air (n mr ) filled pores. All dependences 
of reflectance on wavelength were computed with matrix method [32-35]. Sensitivity of optical sensor is 

—— and depends on porosity P and PS layer thickness L PS , refractive index n ana iyt e (a): 
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where X peak (n analyte ), A peak (n analvte ) - position of maximum of i- peak of PS layer with pores filled with 
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analytes A and B with refractive index n A analyte and n B anal te respectively; 

R( n analyte) ’ R( n analyte) ~ reflection coefficients of structure with pores filled with analyte A, B 
respectively measured at A=const. 
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position and depends on parameters of spectral device that used to measure reflection spectrum: spectral 
accuracy, spectral resolution. 


Accuracy of measurement of selectivity 


A R 


A n 


A-B 

analyte 


depends on reflection coefficient R 


measurements at fixed X and accordingly depends on resolution of analog to digital convertor (ADC) that 
transforms signal from optical transducer into electrical fonn. As fixed X, at which measurements of R for 
set of analytes are performed, rationally to choose the position of maximum of peak A peak (n^ nalyte ) for one 

selected basic analyte from set. Difference R{n^ nal te ) - R{n^ nal te ) for two analytes A and B with slightly 
different refractive indexes can be rather small. 

Noisiness of experimental interference peaks increases with decreasing of integration time of 
collected signal. Noisiness and dissymmetry of peaks complicates detennination of peaks positions and 
amplitude of peak maxima. Procedures of smoothing are used for reduction and elimination of local 
amplitude sharpness and mini peaks. But smoothing can lead to distortion of peak form and displacement 
of peaks positions that can be significant for detection of analytes with slightly different refractive index. 


Colorimetric method 

For finding A c i om (or A com ) and EP the reflectance spectra S(/_) of single PS layer were firstly 
computed using matrix method [32-35]. Tristimulus values X, Y, Z were computed on the base of S(Z) 
using Eq 1-2. Parameters x and y were determined from Eq. 3 to detennine Z dom (or A com ) and EP from 
chromaticity diagram. Spectral dependence of refractive dependence of Si was used for computations [36]. 
Bruggeman effective medium model [37, 38] to link P with the refractive index of PS was used both for 
unfilled pores (filled with air) and pores filled with analyte. 

Dependences of chromaticity coordinates x-y on refractive index of analyte were computed and 
depicted on chromaticity diagram (Fig. 1) for single PS layer with P=50% and Lps=500 mn. Black curve 
ACE (Fig. 1) depicts changing of x and y parameters in accordance with changing of refractive index of 
analyte from n^ lyte = 1 (point A) to n e ^ lyte =1.5 (point E) throw point C(n s ™ yte =1.207). As seen from 



Fig. 1, there are n ana iyt e ranges for which using appropriate reflectance spectra and chromaticity 
coordinates x-y can be calculated and both Ad om and A, C0m can be determined. Both Adom and A com varies in 
wide range of wavelengths and variation depends on initial porosity of layer and L P s. 

To estimate the magnitude of the change of Ad« m , A. CO m, EP in comparison with change of 
interference peak #1 (peak #1 is the first long wavelength interference peak in considered wavelength 
region determined from reflectance spectra) position caused by filling pores analyte with refractive index 
n a naiyte in range 1.31..1.41 corresponding computed dependences are depicted in Fig. 2 for structure with 
P=70% and Lps=500 nm. Shift of position of peak #1 maximum of methanol filled pores comparatively to 
propanol filled pores is (propanol-mcthanol)=Ap Cak |(n propano i)-A peak |(n ine ihanoi) = 14 nm and 

A^fti /An=276. 1 . For analyte pair as water and ethanol computed parameters arc A/Ffj^, =7 nm and 


Al^i /An=244. Shift of chromaticity parameters Adom, A CO m and EP as characteristic of whole spectra are 
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Fig. 2. Computed dependence of shift of dominant , complementary wavelength X'f ' n , the excitation 
purity EP s hi ft , position of diffraction peak on refractive index of analyte n ana iyt e inside the pores. 
Vertical lines depict the value of refractive index for methanol, water, ethanol, propanol. 
























Table 1 


Shift of chromaticity parameters k d0 m, A com and EP for different pairs of analytes 


Parameter 

Pair of analytes 

Analyte A 

Analyte B 

Analyte A 

Analyte B 

methanol 

propanol 

water 

ethanol 


13 

7.6 

ACo.nm 

10.2 

5.2 

A<jL,nm 

6.8 

3.7 

A«, nm 

17.9 

10.5 

A«.nm 

15.7 

8.9 

Ep s hift > nm 

8.8 

5.16 

AAZ/An 

351 

362 

A«/A» 

307.8 

306 

EPshif, 

172.5 

178 


Determined values from Fig. 2 are depicted in Table 1. Comparison of shift of tree interference 
peaks AA s ^ kl 2 3 determined from reflectance spectra calculated by matrix method [32-35] with shift of 

A/l^f , A/tf'f and EP shift shows that parameters k dom , A, com; EP can be used for reflectance spectrum 
characterisation and characterisation of changing of spectrum caused different analyte filling of pores. 
Magnitude of parameters AT^'f , A/tf'f and EP shift that proposed for reflectance spectra characterization is 

larger then magnitude of interference peak shift AA^ kl that traditionally used for characterization of 
spectrum with interference maxima and minima. 

Numerical analysis 

Dependences of kd om , A com and EP on porosity P and refractive index of analyte that fill pores were 
computed. Computed in accordance with Eq. 7 shift of k d om(n ana iyte), A, CO m(n ana iyte), EP(n ana iyt e ) for analyte 
filled pores relatively to values /-dom(n air ), A com (n air ), EP(n a i r ) for unfilled analyte pores (Lps=500nm) are 
depicted in Fig. 3, respectively. There are some ranges of P (for fixed L P s) and n ana iy, e where A dom can not 
be defined due to definition of kd om - These areas are marked with square filling. For PS layers with low 
porosity P<30% A, dom practically isn’t depend on n ana iyt e . Value of increases with P increasing and 

can reach 170 mn for certain n ana iyt e . 

There is limit of maximum value of porosity that can be experimentally fabricated due to 
mechanical strength of PS layers of high porosity caused by internal mechanical tension of porous layer 
through oxidation in air [39]. From practical point of view it is more important to consider difference of 
Adom or A/l s jf for defined set of analytes with defined range n^ alyte ...n™* lyte , where n'Z alyte >n air =l. Vertical 


















lines on Fig. 3 show refractive index range n water = n™™ lyte ...n™Zi yt e = n ethanol that is analyzed later. Values of 
I ^doL I» I K'om I > I EP shifl | increases for higher values of P. 
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Fig. 3. Computed dependence of shift of dominant wavelength 
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^oJ n analyte) ( a )> complementary 


wavelength A com {n amlyte ) (b), EP(n analyte ) (c) on layer porosity and refractive index of analyte n ana iyt e 
inside the pores relatively to A dom (n air ) , A com (n air ) , EP(n ajr ) of air fdled pores, respectively. 


Dependence of Adom, A com . EP on P, n ana iyt e gives 2D matrix as result of computation (Fig. 3). 
Values Adom, Acom, EP also depend on L P s. Dependences of difference 
(ethanol - water ) = X dom (n ethanol ) - A dom (n water ) , AA com (.ethanol - water ) = A com (n ethanol ) - A com (n water ) 
on P for water (n wat er = 1.3236) and ethanol (n e thanoi= 1.3558) were considered to simplify further analysis 
and data visualisation. Factually computed dependences are dependences of differences of profiles at 
fixed n water and n et h an oi specified as line on 3D plot in Fig. 3. 

From practical point of view, for sensor application to receive higher sensitivity AA dom / An , 

A A com / An , for selected range An the optimal parameters P and L P s that detennines maxima |AA ( i„ m |, 
|A^com|, should be found (Fig. 4). 









































































Fig. 4. a) Computed dependence of local maximums of difference A dom (n elhanol ) - X dom (n water ) on porosity 
P and PS layer thickness L PS ; b) computed dependence of shift of dominant wavelength X dom (n aimlye ) and 
excitation purity EP(n ana iyt e ) for analyte fdled pores in reference to X dom (n air ) and EP(n air ) of air fdled 

pores respectively on refractive index of analyte n ana iyt e inside the pores for PS thickness L PS =800 nm and 
P=69.5% (point A, Fig. 4a). 

Local maxima represented on Fig. 4a with circles, local minima with squares. For sample with 
certain L P s there can be several local maxima and minima at discrete values of P. Analysis of dependence 
of |AXdom| on discrete values P and L P s shows that position of local maxima |A/.d,, m is described with 
uninterrupted curves on plane of P and L P s variables. These continuous curves that depicts dependence of 
maxima of | X dom (n ethanol ) - X dom (n water ) | on L PS and P are displayed as dotted lines in Fig. 4a. So, to receive 

maximum response of sensor as maximum values of X dom (n ethanoI ) - X dom (n water ) for analyte as ethanol 
and water one should choose point on these curves for desired L P s and P. Analyze of X dom (n ethano[ ) - 
X dom ( n W ate r ) values along dotted lines on Fig. 4a shows presence of maximums: 52.5 nm (L PS =400 nm, 

P=53%), 37 nm (L P s=800 nm, P=56.5%), 102.7 nm (L^^OO nm, P=69.5%), 32 nm (L P s=900 nm, 
P=78%)... 

Computed dependence of maximum shift of dominant wavelength A dom (n anal ) and excitation 
purity EP(n ana iyte) in reference to X dom (n air ) and EP(n air ) respectively on refractive index of analyte n ana iyt e 
for structure with L P s=800 nm and P=69.5% (point A, Fig. 4a) is shown on Fig. 4b. For pair water and 
ethanol I X dom (n ethanol ) - X dom (n water ) | is 102.7 nm and EP(n ethmol ) - EP(n water ) =1.21%. For pair propanol 
and ethanol \A dom (n propano/ )-A dom (n ethanol )\ is 4.8 nm and EP(n pwpanol )~ EP(n ethanol ) is 0.73%. These values 

for ethanol-propanol are less comparatively with values for ethanol-water because parameters L P s and P 
are optimal for analyte pair as ethanol-water. Optimal parameters P for any pair of analyte as 








































analyte A-analyte B can be detennined from Fig. 3 for structure with L P s=500 nm. To receive optimal 
parameters L PS dependences as depicted on Fig. 3 should be computed for different values of L PS . 

Adsorption/desorption kinetics 

Experimental dependence of reflection coefficient on time for PS sample (L P s=600 nm, P=50%) 
with pores filled with ethanol depicted in Fig. 5. Initial time to is moment of beginning of deposition of 
ethanol on the surface of PS sample. There are transition processes related with analyte deposition and 
pores filling during time interval to..fi. Time range ti < t < t 2 corresponds to the state when pores are filled 
with analyte and analyte is evaporated from surface. Positions of peaks maximums practically doesn't 
change during analyte evaporation from sample surface due pores are still filled with analyte. Thickness 
of liquid analyte layer on the surface of sample decreases with time due to analyte evaporation and this 
results in jitter of widths of interference peaks maximums in time range fi < t < t 2 . The invariance of the 
position of the interference peaks and the jitter of the peak width upon evaporation of the analyte from the 
PS surface was shown as a result of computer simulation in work [35, 40], Desorption of analyte 
molecules from pores begins at t 2 and is described by time range t 2 ..t 3 . Region t>t 3 accords to state when 
PS layer reflectance spectrum came back to initial state, i.e. to full analyte desorption both from PS 
surface and pores. Reflectance spectrum has non monotonic behavior during desorption in region t 2 ..t 3 . 
Reflectance spectrum shape change at t>t 3 comparatively to initial spectrum at t<to can be caused by 
presence of un-desorbed molecules of analyte inside pores and depends on interactions of analyte 
molecules with molecules that form pores. 

Filling of pores with liquid analyte or in other words replacement of air inside the pores with an 
analyte results to clear shift of interference peaks due to large difference of refractive index of air and 
liquids. But difference of interference peaks position can be small for analytes that have small refractive 
index difference. Detectable interference peaks position difference depends on spectrometer parameters, 
noise level, parameters of porous structure. 

Reflectance spectrums for PS structure (L PS =500 nm, P=50%) with pores completely filled with 
water and ethanol after completion of all transient processes are shown in Fig. 6. Integration time of 
measured separate spectrums was 50 ms to measure kinetics of adsorption/desorption processes. 
Experimental spectrums are noised and were smoothed to detennine positions of peaks maximums. Shifts 
of positions of peaks maximums in observed wavelength range are A shifi =9 nm (long wavelength peak), 

?i shif =2.8 nm, fi shift =1.5 nm (short wavelength peak). Peak shift A shift is larger for peaks in long 

wavelength part of spectrum but simultaneous widening of peaks and noisiness complicates peak position 
detennination. At moment t 2 desorption processes from pores starts. These processes accompanied by the 
reallocation of analyte inside pores and on the surface of sample that result to complicate form of 
reflectance spectrum, splitting of peaks (Fig. 6 inset). Presence of complex shape peaks with several local 
maxima makes impossible to determine the exact position of the peaks especially if analytes with slightly 
different refractive indexes are considered. 
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Fig. 5. Experimental dependence of reflection coefficient on time for sample with one PS layer (Lps=600 

nm, P=50%) filled with ethanol. 



Fig. 6. Experimental dependence of reflection coefficient on wavelength for sample with PS layer 
(Lps=500 nm, P=50%) filled with water and ethanol in stationary state and transient state (inset). 






























x coordinate 

Fig.7. 1931 chromaticity diagram and shift of experimental sample color chromaticity coordinates x-y 
after pores filling with ethanol (L PS =500 nm, P=50%). 

Chromaticity coordinates x-y were computed on the base of dependence of reflectance spectra 
( Xnm =420 nm to A max =780 nm) on time for PS sample (L P s=500 nm, P=50%) with pores filled with 
ethanol as analyte (Fig. 7). Wavelength X mm =420 nm was chosen as lower limit for measurement to avoid 
spectra noiselessness at a< 420 nm during analysis of kinetics of adsorption/desorption processes. 
Increased noise is observed in short wavelength range due to decreasing of the grating efficiency of used 
CCD spectrometer and small integration time required for measurement of the quickly changing spectrum. 
Point A represents reflectance spectrum for pores filled with air. Point B represents reflectance spectrum 
for pores filled with ethanol. For unfilled pores Adom(n ai ,-)=563 nm and for pores filled with ethanol 
A-dom(n e thanoi) = 574 nm. All chromaticity coordinates and chromaticity parameters /Com, A. com , EP were 
computed for measured reflectance spectrums limited with range /. mm =450 nm to /- max =780 nm to avoid 
the influence of noisy low wavelength range peaks. 

Kinetic of adsorption/desorption processes for ethanol and water expressed as dependence of 

Alom ( ^analyte ) ^ dom ( ^ analyte ) ^ dom (O and \ om (n analyte ) on time is depicted on Fig. 8a for PS sample 

with L P s=500 nm and P=50%. Dependences for ethanol and water were aligned relative to the time to of 
analyte deposition. There are transition processes related with analyte deposition and pores filling in the 
range to-.ti- Clear shift of A<iom and X com is observed after completion of transition processes and complete 
filling of pores with analyte at moment fi. Shift '^,( n et h ano i ) °f ^dom(n e thanoi) *' or sample with pores filled 
with ethanol relatively of kdom(n a ir) for unfilled pores is 10.3 nm. For pores filled with water 








%do!n Ohvter ) = KoM.cer) ~ K om ("air ) =8-7 nm. Difference of X™ for ethanol and water is 
A 4, m (ethanol- water) = &n ethanol )-A s d Z(n waler ) = 1.6 nm and ratios are 


A^dom(ethanol- water)/(n et hanoi _ n W ater) = 43 and AA com (ethanol-watcr)/(n el hanorn W at cr )=446.2. Dependences of 
AW=WC 2 H 5 0H-H 2 0 mixture)-7, dom (H 2 0) and AX CO m=7.com(C2H 5 OH-H20 rnixture)-A. com (H 2 0) on 
ethanol/water volume ratio can help to distinguish mixtures of water and ethanol (Fig. 8b). 




Fig. 8. a) experimental dependence of shift of dominant wavelengths A s d *(n ethanol ) , (n m , ter ) and 
dependence of complementary wavelengths X com (n ethanoI ) , A com {n water ) on time after deposition of ethanol 
and water on sample surface and pores fdling (L P s=500 nm, P=50%); b) experimental dependence of 
A^dom, A/, com on ethanol/water volume ratio (Lps=500 nm, P=50%). 

As was shown before whole spectrum can be characterized with ^domor A com . Value of AA, com shifts 
up to 442 nm after filling of pores with water and shifts up to 457 nm after filling of pores with ethanol. 
Difference of A com for pores filled with ethanol and water is 

AA com (ethanol - water) = k com (n ethanol ) - \ :om (n water ) = 15 nm. Value A com isn’t defined for time t<t 0 and for 
times after desorption of analytes from pores starts due definition of /. com . Ratio 
A^com(ethanol-water)/(n eth anoi _ n wa ter) = 403. Desorption of analytes begins at ti; and t" and ends at t^and 
t" for ethanol and water respectively. Different time of analyte desorption from pores At* s anol < Af“ crl is 

explained with the difference in wettability and interaction of water and alcohol molecules with silicon 
walls in the pores 

Kinetic of adsorption/desorption processes for ethanol, water and solution ethanol-water (1:1) was 
investigated and depicted as dependence of ^(n atmb J=A dom (n ana i yle ) - A dom (n air ) on time on Fig. 9 for PS 
sample with Lps=600 nm and P=50%. Dependences for all analytes were aligned relative to the time to of 
analyte deposition start. After completion of transients processes at fi value A d ^ n stays practically constant 
till beginning of desorption processes at t 2 . For investigated analytes there is 




























relation X^(n water ) < ^!Un ethanol - water u) < ^(n ethanol ) ■ This relationship is confirmed by computer 
simulation on Fig. 4b were l? d f‘(n water )< ^(n ethanol ). 



Fig. 9. Experimental dependence of shift of dominant wavelengths X^{n ethanol ) , A s d h *(n water ) , 
KonX n ethanol-water) 011 l ' mc after deposition of ethanol, water and ethanol-water solution (1:1) on sample 
surface and pores filling (L P s=600 nm, P=50%). 

Conclusion 

CIE 1931 color space based parameters as k dom , A com , EP can describe kinetics of 
adsorption/desorption processes of analyte in PS single layer structure. Colorimetric analysis based on 
parameters as A, dom or A, com characterizes spectrum changing as whole contrary to determination of 
changing of parameters of separate spectrum features as interference peaks. Computed dependences of 
kdom, k c „m, EP on porosity (P), PS layer thickness (L P s) shown the presence of sets [P, L P s] where 
maximum selectivity is observed for analytes with refractive indexes in the range n m i n ...n max . Time 
dependences of k dom , A com show complex behaviour at adsorption/desorption of ethanol, water and 
ethanol-water solution. Ratios Ak dom ( ethanol-water)/(n et hanoi _ n wat er) = 43 and 

A/. com (cthanol-watcr)/(n c th an onn walcI )=403 were experimentally obtained for ethanol and water. Using 
parameters as /. dom , A, com , EP is attractive approach to reduce data matrix dimension or number of sensor 
elements in sensor matrix of electronic tongue/nose based on single PS layer. 
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